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Confinement in non-Abelian gauge theories, such as QCD, is often explained using an analogy 
to type II superconductivity. In this analogy the existence of the "Meissner" effect for quarks 
with respect to the QCD vacuum is an important element. Here we show that using the ideas of 
Abelian projection it is possible to arrive at an effective London equation from a non-Abelian gauge 
theory. (London's equation gave a phenomenological description of the Meissner effect prior to the 
^SJ ' Ginzburg-Landau or BCS theory of superconductors). The Abelian projected gauge field acts as 

' the E&M field in normal superconductivity, while the remaining non-Abelian components form a 

, gluon condensate which is described via an effective scalar field. This effective scalar field plays a 

04 ' role similar to the scalar field in Ginzburg-Landau theory. 

P I PACS numbers: 12.38.Aw, 12.38.Lg 
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I. INTRODUCTION 

> 

■ One of the difRcult aspects of QCD is that it is a nonlinear theory. At the classical level one writes down the field 
[ equations of the nonlinear theory, and by inspiration or numerically solves the field equations. For non-Abelian gauge 

■ theories this leads to interesting classical field configurations such as monopoles, merons, instantons etc. In the quan- 
' tized version of these theories there are suggestions of interesting configurations like those listed above. For example, 
, the chromoelectric flux tubes that are hypothesized to exist between quarks in the dual, color superconductivity 

picture of confinement. 

Quantizing nonlinear field theories presents additional challenges and questions. In the case of general relativity 
no complete quantum version of the theory exists. In the case of QCD a quantum version of the theory exists, but 
O i' because the QCD coupling constant is large, the standard perturbative quantum field theory techniques do not apply 
^ . in the low energy regime. This means that getting numerical predictions out of QCD is difficult. At present there is 
no universal analytical technique for handling detailed questions about QCD (this is in contrast with QED and the 
theory of the electroweak theory where the techniques of perturbative quantum field theory supply a general method 
for calculating results analytically). The nonperturbative aspects of QCD are dealt with using numerical techniques 
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' In the I950's Heisenberg and co-workers proposed a scheme for investigating quantum, nonlinear spinor fields 
which yielded interesting results concerning the spin and charge of the electron. This nonperturbative quantiza- 
tion technique distinguished between field operators of interacting fields and field operators of noninteracting fields 
(Heisenberg's ideas were developed before the widespread adoption of the path integral techniques for quantizing 
field theories and as such are given in terms of field operators). The algebra obeyed by the noninteracting field 
operators is just the standard canonical commutation relationships (see for example pg. 86 of ref. |^). The algebra 
for the interacting fields, on the other hand, is determined from the Green's functions for the field operators of the 
nonlinear fields, G{xi,X2 ■ ■ . ,Xn) = {Q\A{xi), A{x2), ■ . . A(x„)|Q), where A{xi) is the field operator at position Xi; 
G{xi, X2 ■ ■ ■ , Xn) is the Green's function; \Q) is a quantum state. Conversely if the G{xi, X2 ■ ■ . , Xn)s are known then 
this gives the quantum states \Q). The Green's functions are determined from the field equations in the following way: 
one starts with the classical field equations and turns the classical fields into field operators. One uses these operator 
field equations to determine the lowest order Green's functions. Next one applies various order derivative operators 
on these operator field equations to obtain equations which have higher order powers of the field operators, and which 
yield the higher order Green's functions. In this way an infinite set of coupled, differential equations are constructed 
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which connect all orders of Green's functions. In general it is not possible to solve this system of differential equations 
analytically, so approximation techniques must be used. 

Here we apply a version of Heisenberg's method to non-Abelian gauge theory. In our treatment we assume that 
the non-Abelian gauge fields can be separated into two classes: stochastic, disordered fields and ordered fields. This 
is similar to Abelian Projection ideas where the diagonal gauge fields associated with the Abelian subgroup are 
ordered (an example of this ordering would be a non-zero expectation, {Af''"-^) ^ 0, inside a flux tube stretched 
between quarks) and with the oflF-diagonal gauge fields being disordered so that (A'^f^) = outside of the flux tube 
(but {{A'^f^)'^) 7^ outside the flux tube). This picture also has some common points with the work of Nielsen, Olesen 
and others Q - |^ on the response of the QCD vacuum to a homogeneous color magnetic field, H. In ref. Q it was 
shown that a homogeneous SU(2) magnetic field lowered the energy of the vacuum, but was unstable. This instability 
was removed by the formation of domains These domains took the form of SU(2) magnetic flux tubes which 
formed a random SU(2) magnetic quantum liquid ||] with the property that (H) = and {H"^) ^ 0. The difference 
with our approach outlined above is that we assume the expectation values are in terms of the gauge potentials rather 
than the magnetic fleld - (A'^^^) = and {(A'^f^y) ^ 0. Also, in line with the ideas of Abelian Projection, the 
diagonal, Abelian potentials (Af*"^) and the off-diagonal, non-Abelian potentials (A'^^f ) play different roles in the 
present paper. 

In applying these ideas to an SU(2) gauge theory we assume that the off-diagonal components form a condensate 
which can be described by an effective scalar field similar to that in Ginzburg-Landau theory |^ . This is also related 
to work by Cornwall ||l^ where a connection is established between a scalar field and the expectation of the Yang- 
Mills field strength tensor squared - (j){x) ^ (TrGf^^). In the present work we follow the Abelian Projection ideas 
and assign different roles to the Abelian and the off diagonal, non-Abelian gauge potentials. Only the off-diagonal 
gauge potentials are assumed be involved in the condensate. With respect to this condensate the field equation for 
the Abelian component then takes the form of the London equation, so that the Abelian field develops a mass as it 
penetrates into regions characterized by the a non-zero value for off-diagonal condensate. The mass of the Abelian 
field and the assumed effective mass of the condensate are not the same, which is in contrast to refs. jll] where 
all the SU(2) gauge potentials play an equivalent role, leading to a single common effective mass for all the gluons. 



II. SEPARATION OF COMPONENTS 

In this section we follow the conventions of Ref. Starting with the SU(N) gauge group with generators we 
define the SU(N) gauge fields, A,, = A^T^ . Let G be a subgroup of SU(N) and SU{N)/G is a coset. Then the gauge 
field A^ can be decomposed as 

A^ = A^T^ - a^T"^ + A^^'T", (1) 
e G and A™ e SU{N)/G (2) 

where the indices a,b,c. . . belongs to the subgroup G and m, n, . . . to the coset SU(N) /G; B are SU(N) indices. Based 
on this the field strength can be decomposed as 

j'l^.T^ ^t;:,t^ + t;it"' (3) 

where 

= + e G, (4) 

= d^al - d^al + r'^alal G G, (5) 

= r^^^A^A: eG, (6) 

= + C;", e SU{N)/G, (7) 

Fjr, = d^A"J - d^A"; + r'^^'AlAl G SUiN)/G, (8) 

C-m ^ frnnb ^^n^b _ ^n^b) ^ SUiN)/G (9) 

where /^^*^ are the structural constants of SU(N). The SU(N) Yang-Mills field equations can be decomposed as 

{F""^" + G'"'"') = -/""'' [A'J + ^^f"") - al (i^"^'' + G"'"')] (11) 

where • • ]" = d^[- • • ]° -I- /"'"^a(^[- ■ ■]'^ is the covariant derivative on the subgroup G and D^[- • • ]™ = d^[- ■ ■ -I- 
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Specializing to the SU(2) case we let SU{N) SU{2) , G U{1), and /-^^^ -> e^^^ . 
a = 3 and consequently m,n = 1,2, our classical equations become 

Since G = [/(I) we have = d^. 



Setting the indices as 

(12) 
(13) 



III. HEISENBERG QUANTIZATION 



M and 



fields by field operators 



A 



In quantizing the classical system given in Eqs. (jlj) - ( |13| ) via Heisenberg's method one first replaces the classical 

™. This yields the following differential equations for the operators 



^3mn 



a:: 



G'^ 



(15) 



These nonlinear equations for the field operators of the nonlinear quantum fields can be used to determine expectation 



values for the field operators and A™ (e.g. 



where (• 



(01 • • • \Q) and \Q) is some quantum state). One 



can also use these equations to determine the expectation values of operators that are built up from the fundamental 
operators and A™. For example, the "electric" field operator, = c^oOz — S^ao giving the expectation (i?z). The 
simple gauge field expectation values, (^^(a;)), are obtained by average Eqs. ( p^ ( psj ) over some quantum state \Q) 



Q 



3mn 



G" 



0) = 0, 



(f™^'' + G 



+ € 



3mn 



A] 







(16) 



(17) 



One problem in using these equations to obtain expectation values like (A™), is that these equations involve not only 
powers or derivatives of {A^) {i.e. terms like or dadp{A^)) but also contain terms like C/™ = (A™ A"). 

Starting with Eqs. (p^)- ( p7| ) one can generate an operator differential equation for the product 4™ A" thus allowing 
the determination of the Green's function . However this equation will in turn contain other, higher order Green's 
functions. Repeating these steps leads to an infinite set of equations connecting Green's functions of ever increasing 
order. This construction, leading to an infinite set of coupled, differential equations, does not have an exact, analytical 
solution and so must be handled using some approximation. 

Operators are only well determined if there is a Hilbcrt space of quantum states. Thus we need to ask about the 
definition of the quantum states \Q) in the above construction. The resolution to this problem is as follows: There is 
an one-to-one correspondence between a given quantum state \Q) and the infinite set of quantum expectation values 
over any product of field operators, Q™J!^y {xi, X2 . . . ) = {Q\A^^{xi)A^{x2) . ■ . \Q). So if all the Green's functions - 
t/JJ™.'.' (xi,X2 • . . ) - are known then the quantum states, \Q) are known, i.e. the action of \Q) on any product of field 

operators A™ (a; i) A" (0:2) ... is known. The Green's functions are determined from the above, infinite set of equations 
(following Heisenberg's idea). 

Another problem associated with products of field operators like A™{x)A2{x) which occur in Eq. (p^ ) is that the 
two operators occur at the same point. For non-interacting field it is well known that such products have a singularity. 
In this paper we are considering interacting fields so it is not known if a singularity would arise for such products 
of operators evaluated at the same point. Physically it is hypothesized that there are situations in interacting field 
theories where these singularities do not occur {e.g. for flux tubes in Abelian or non-Abelian theory quantities like 
the "electric" field inside the tube, (F") < 00, and energy density s{x) = ((F^)^) < 00 are nonsingular) . Here we 
take as an assumption that such singularities do not occur. 

We now enumerate our basic assumptions: 

1. After quantization the components A™(a;) become stochastic. In mathematical terms we write this assumption 
as 



{Ai:{x)) = o 



and 



{ai:{x)a:{x)) = -^(x)^™"77^. 



(18) 



where (p{x) is some scalar field, rj = {+1, —1, —1, —1}. This would give a problem with the time components in 
that (A™ A™) < 0. Thus to deal with this we also assume that the fields are static and have no time component. 



I.e. 
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0. 
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2. The components a° of the subgroup G can have some order so that certain expectation values can have non-zero 
values, for example 

(H^) = ((V X a),) ^ 0. (19) 

Such conditions are meant to imply that (or certain quantities derived from it) develops a non-zero expectation 
value for some non-trivial, non- vacuum boundary conditions {e.g. the presence of external quarks). Such 
conditions are not connected with vacuum states since this would imply a violation of the Lorentz symmetry of 
the QCD vacuum. 

3. The gauge potentials a'^ and A™ are not correlated. Mathematically this means that 

(/«).9(^™)) - (/«)) (giA-;)) (20) 

where /, g are any functions 

These assumptions are a variation of the Abelian Projection ideas, since there the SU(N)/G components of the gauge 
fields are suppressed. The characterization of the off-diagonal fields as stochastic is a result of the first part of Eq. (|lq), 
{A]^{x)) = 0. The second part of Eq. is related to some recent work which demonstrates the physical 

importance of the expectation value of the square of the non- Abelian gauge potential to the dynamics of non- Abelian 
field theory. The surprising thing about this is that the non- Abelian gauge potential (and its square) is gauge variant, 
and one would think that physical quantities should only be constructed from gauge invariant quantities. In previous 
work Q ||T^ one had conditions similar to the first assumption above, but in terms of the expectation values of 
the Yang-Mills field strength tensor and its square - (G^jy) — and (G^^) ^ 0. One way of looking at the condition, 
{x) (x)') = —Lp{x)S""^ri^^, is that it represents the condensation of the off-diagonal SU(2) gluons into effective 
scalar fields, ip{x). This provides a physical motivation for a connection of the present work to the Ginzburg-Landau 
model of superconductivity. In Ginzburg-Landau theory the scalar field represents a condensation of electrons i.e. 
the Cooper pairs. This association between the expectation of the square of the off-diagonal gauge potentials with a 
scalar field is also similar to ref. ||l^ except there the association was between (Gj;^) ^ and the scalar field. 



IV. LONDON'S EQUATION 



In this section we want to show how London's equation emerges from Eqs. (14)-([15[) under the setup outlined 
above. London's equation describes the Meissner effect in ordinary superconductivity. Showing that the same equation 
emerges from a quantized non-Abelian gauge theory gives support to the dual superconducting picture of the QCD 
vacuum. Because of the stochastic assumption above we will not be interested in the off-diagonal components of the 
gauge fields, (A™). Thus we will not worry about Eq. (^5|) which is the equation that determines these off-diagonal 



components. The Abelian field Om is determined from Eq. (14) which is linear in a^. Because of this we take the 



Abelian gauge field as classical ||lj| . This leads to the following equation 



9.1 



((A^J'F"^"^) + {A1'G''^"')) . 



(21) 



Note the Abelian term, 0'^", is treated classically while the remaining terms which involve combinations of the off- 
diagonal fields are treated as quantum degrees of freedom via the expectation values. To calculate these expectation 
values we take, as a first approximation, the scalar function of Eq. (nS) as a constant, i.e. (p{x) = (po 



{A^ix)A:ix)) = -^oS""'r,,, 



Then this gives 



The next term is 



(v4™G"^'') 



(a^a:) = 0, 



((A™v4f^) a" - {A'^AP") a'') = -3(/:>oe^""'a''- 



(^^rnpni^u^ = (^1^9''^"'') - + e^P" {A"^ AP^" A"") 



(22) 

(23) 
(24) 

(25) 



For the disordered, non-diagonal components we will set (^"^^ {x)A™^ {x) . . . {x}) = if n is odd. For the other 
terms in the right-hand-side of Eq. ( p5| ) we note that 

9„ (A-; a:) = (9„A';M'J) + (A^O^A:) = O, (26) 
(d^A^^A:) = - (A^'d^A:) . 
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For these stochastic, non-diagonal components the last expression should not depend on the order of the indices (m, n) 
and i.e. {daA]^A'^) = (d^A'^A]^) = Using this with Eq. ^ gives 



(d^A^A:) = (A^do^A:) = 0. (27) 
Putting all this together gives from Eq. ( pl| ) 

a^c/.A"' = Qipoa". (28) 

applying the Lorentz gauge condition, di^a'^ — 0, then yields 

d^.d^a" = -Gifioa". (29) 

This is London's equation for the U(l) ordered phase in the presence of disordered SU(2)/U(1) phase. To demonstrate 
how this leads to a Meissner-like effect for the U(l) gauge field, a^, we take half of all space as being filled by the 
stochastic phase {e.g. (f{x) ^ ipo ^ for y > and (p{x) — for y < 0). In this case the Abelian gauge field has only 
a dependence on y, a^{y), and Eq. (p^ ) becomes 



&Lpoa^ (30) 



dy"^ 

which has the solution 

= ao^^e-^y. (31) 

Thus the magnetic field Hz = ffoze"^^'^"^ is exponentially damped as it penetrates the region with the stochastic 
phase. 

From eqs. (g)-(|l|) the effective mass of the Abelian field is rrie f f — \/Q(po- On the other hand eq. ( p2| ) (up to 



a group factor of 2/3) is similar to the relationship given in ref. |11| (see eq. (32) of that reference) between the 
effective gluon mass and the expectation of the square of all the gauge potentials. Thus from eq. ( p^ ) we find that in 
our model the effective mass of the SU(2) gluons associated with the off-diagonal gauge potentials is different from 
the effective mass of the gauge boson associated with the Abelian gauge potential. The difference in masses between 
the condensate represented by the scalar field and the gauge boson is also found in spontaneous symmetry breaking 
of a gauge symmetry. For example, consider the Ginzburg-Landau Lagrangian with an Abelian gauge field, A^, and 
a complex scalar field, if, describing the condensate 

Cgl = {D^ip){D^ipr - m^\cp\^ - X\ip\^ (32) 



where -D^ = d^^ + ieA^. The condensate has a mass of m while the gauge boson A^ will acquire a mass of w 

In our case the condensate comes from the same set of SU(2) gauge fields as the Abelian gauge field. The different 
behavior/roles of the Abelian and off-diagonal, non- Abelian gauge fields results from using ideas similar to Abelian 
Projection through our first assumption given in eq. (^8|) above. 



V. CONCLUSIONS 



In this paper we have shown how the London equation emerges from a non- Abelian gauge theory by combining ideas 
of a nonperturbative quantization technique pioneered by Heisenberg and co-workers, with ideas similar to Abelian 
Projection. The importance of this is that the London equation gives a phenomenological description of the Meissner 
effect in superconductors, and the vacuum of some non- Abelian gauge theories {e.g. QCD) is often modeled as a dual 
superconductor in order to explain confinement. In our approach we split the gauge group (SU(2) in our case) into a 
subgroup (U(l) in our case) and the coset space (SU(2)/U(1) in our case). The gauge bosons associated with the coset 
SU(2)/U(1) were taken to be in an disordered, stochastic phase, (A™(x)) = 0. Mathematically this statement was 
contained in Eq. (|l^) where the scalar field can be compared to the scalar field in the Ginzburg-Landau treatment of 
superconductivity. In the Ginzburg-Landau model the scalar field represents a condensation of electrons into Cooper 
pairs. In our work the scalar field can be thought of as a condensation of gluons. Just as the E&M field is excluded 
from the superconductor, so in our example the diagonal, Abelian gauge field is excluded from the disordered phase. 

There is a difference between Abelian Projection and the treatment in the present paper. In Abelian Projection 
the off-diagonal components are constructed by applying gauge fixing, but in our case they emerge from applying the 
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three assumptions given in section III to the dynamical equations. In the first approximation we have neglected the 
dynamical behavior of the stochastic phase, by setting Lp{x) = <^o- In this way we obtained an equation for the Abelian 
components of the gauge field which was similar to London's equation for the vector potential in superconductivity 
theory. Higher order approximations in the above procedure would result in higher order powers and derivatives of 
Lp. This would hopefully lead to dynamical equations for Lp{x) similar to the field equations which result from the 
Ginzburg-Landau Lagrangian given in eq. (^2|). (Note in this regard that there are two scalar fields in eq. (|3^), since 
there tp is complex, and in eq. ( p^ ) there are also effectively two scalar fields coming from m = n = 1 and m = n = 2). 
In this case one would be able to construct Nielsen-Olesen flux tube solutions |lq| , which would be very suggestive 
toward making a firm connection with the dual superconducting model of QCD. Such a construction of an effective 
Ginzburg-Landau equation for ip would be important in bolstering the claim of a connection between our approach 
and the dual superconducting model of the QCD vacuum. 
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